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Abstract 



Modifications in the spontaneous emission rate of an excited atom that are caused by extinction effects in a nearby dielectric medium 
'are analyzed in a quantummechanical model, in which the medium consists of spherical scatterers with absorptive properties. Use of 
^—1 the dyadic Green function of the electromagnetic field near a a dielectric sphere leads to an expression for the change in the emission 
rate as a series of multipole contributions for which analytical formulas are obtained. The results for the modified emission rate as 
a function of the distance between the excited atom and the dielectric medium show the influence of both absorption and scattering 
C\ processes. 
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' The emission rate of an excited atom is modified if the elec- 
' P_| ' tromagnetic properties of its surroundings differ from that of 
Q^^yacuum For an atom in front of a dielectric medium fill- 
^ I ^ ing a half-space, the rate varies with the distance between the 
^ atom and the medium |01 - iSl- Usually, the medium is taken 
to be homogeneous on the scale of the atomic wavelength, so 
^^ that its electromagnetic properties are fully described by a sus- 
I I ceptibility, which does not vary appreciably on the scale of the 
wavelength. In general, it will be complex to account for ab- 
sorption and dispersion. For such a configuration modifications 
of the atomic radiative properties have been confirmed experi- 
mentally a few years ago [9]. 

If the structure of the medium cannot be neglected, seat- 
ed tering effects play a role as well, so that extinction in such a 
medium is driven by both absorption and scattering. Extinction 
by scattering in material media is quite common, owing to the 
presence of impurities and defects. The interplay between the 
\ [ two types of extinction in atomic decay rates can be investigated 
J> in a model in which both of these features occur simultane- 
ously. In a recent paper 1 10] we have studied the change in the 
decay rate of an atom in the presence of a medium consisting 
d of non-overlapping spheres that are made of absorptive dielec- 
tric material. The spheres are distributed randomly in a half- 
space, with a uniform average density. In order to describe the 
absorptive dielectric material of the spheres in a quantumme- 
chanically consistent way a damped-polariton model has been 
employed |11]. By introducing an effective susceptibility for 
the composite medium and after a detailed analysis of surface 
contributions, we could derive an asymptotic expression for the 
change in the emission rate at relatively large distances between 
the atom and the medium. 

In the present paper we take a somewhat different approach 
so as to derive an analytic expression for the emission rate that 
is valid for all distances between the atom and the medium. We 
shall start from exact expressions for the electromagnetic Green 
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function in the presence of a dielectric sphere of arbitrary ra- 
dius. As before, the absorptive dielectric material of the spheres 
will be described by a damped-polariton model. For simplicity 
we shall assume that the density of the spherical scatterers in 
the medium is low, so that multiple-scattering effects can be 
neglected. 

2. Spontaneous emission in the presence of absorbing di- 
electrics 

In the damped-polariton model an absorptive linear dielec- 
tric medium is described by a polarization density that is cou- 
pled to a bath of harmonic oscillators with a continuous range 
of frequencies [11]. The Hamiltonian of the damped-polariton 
model can be diagonalized exactly, as has been shown both for 
the case of a uniform dielectric [ 1 1 ] and for a dielectric with 
arbitrary inhomogeneities lll2ll . Diagonalization of the Hamil- 
tonian in the general non-homogeneous case yields 



/-I 



dtj Hcu C ' (r, w) ■ C(r, w) , 



(2.1) 



with annihilation operators C(r, Lui) and associated creation op- 
erators. The electric field can be expressed in terms of these 
operators as lll2ll : 



E(r) 



d(AjiE(r,r' ,Lo) ■ C(r',(y) + h.c. , 



(2.2) 



with a tensorial coefficient: 

up' Iti Ime(r', u + iO) 



f£(r,r',w) = -/-^ 
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G(r,r',w-H/0) .(2.3) 



Here s is the complex local (relative) dielectric constant, which 
follows from the parameters of the model. Furthermore, G 
is the tensorial Green function, which satisfies the differential 
equation 



-Vx [VxG(r,r',w + /0)] 

^2 

H--- e(r, a> + iO) G(r, r', co + iO) = 1 5(r ■ 



r'), (2.4) 
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with I the unit tensor. 

The atomic decay rate in the presence of an absorbing di- 
electric follows from the inhomogeneous damped-polariton mod- 
el in its diagonalized form by employing perturbation theory in 
leading order ifioll . It can be expressed as an integral over a 
product of the coefficients (12.3b and suitable atomic matrix ele- 
ments: 



In 



J jdr jdr' j dr"{e\Jair')\g)-iE(r',r,co,y 
■fl(r",r,oj,)-{g\Ur")\e}, 



(2.5) 



with e and g denoting the excited and the ground state of the 
atom, <jL>a the atomic transition frequency, and the tilde denoting 
the tensor transpose. Furthermore, Jfl(r) is the atomic local cur- 
rent density -^e 2i{P(7'w> ^(r-r,)}, with r,- , p, the positions and 
momenta of the electrons and curly brackets denoting the anti- 
commutator In the electric -dipole approximation the atomic 
decay rate can be expressed in terms of the Green function as: 
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^<e|//|g>-ImG(r <^a + iO) ■ {g\fi\e) , 



(2.6) 



with r,, the atomic position and - -e 2,(r, - To) the atomic 
electric dipole moment. The above expression for the decay 
rate of an excited atom in the presence of an inhomogeneous 
absorptive dielectric can be obtained as well by invoking the 
fluctuation-dissipation theorem Iil3nl4i1 . 

3. Green functions 

The Green function in vacuum fulfils the differential equa- 
tion (I2.4l l with e = 1 . It follows from the scalar Green function 
GAx,r',(jj) = exp(/w|r - r'|/c)/(47r|r - r'|) as Go(r,r',w) = 
-[I + (c^/w^)VV]G.5(r, r', a)). Its explicit form in spherical co- 
ordinates is obtained from the expansion of G , in spherical har- 
monics and spherical Bessel functions. The ensuing form for 
the vacuum Green function is lll5[ll6ll 



tX) c 



Go(r,r',../0) = -/.2Z^ 

X [B{r - r') [M(^/r)M^,_„(r') + N<^,(r)Nf,_„,(r') 
+ e(r' - r) [Mf,^(r)M<*l,„(r') + No„(r)N<^'l,„(r')]} 



+k ^ trtr S(r - r') , 



(3.1) 



with k = oj/c the wavenumber, e,- a unit vector in the direction 
of r and 6(r) a step function that equals 1 for positive and for 
negative argument. The vector harmonics are defined as 



Mf,,„(r) = V A [ri//e,,n(r)] , 
N^,„,(r) = riVA[VA[riA^,„(r)]]. 



(3.2) 
(3.3) 



where \l/( „,{r) stands for jt{kr) ¥(^^(0, (p), with jV spherical Bessel 
functions and ,„ spherical harmonics. The superscripts {h) 
in (13. Il l denote the analogous vector harmonics with spherical 
Hankel functions h^P instead of ji. 

The expression ( 13.1b may be checked by substitution in the 
differential equation (12.4b . Differentiation of the step functions 



yields singular terms, which together with the last term lead to 
the right-hand side of (12.4b . 

The Green function in the presence of a dielectric sphere is 
the sum of the vacuum Green function and a correction term. 



For a sphere centered at the origin the latter has the form QIS |- 



c=\ m=-e 
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[B^N<^i(r)N^!jr') + fi?-M<:i(r)M(^l„(r')l 



(3.4) 



for r and r' both outside the sphere. The electric and magnetic 
multipole amplitudes read iflsll - lEoll : 



1 2^+1 < 
( ' ({(+ 1) If, 



(3.5) 



with p - e, m. The numerators and denominators are given as 
N'e=Efe{q)ie{q')- jt{q)fi{q'). 



NT =fe(q)je(q')-jc(q)Mq'), 
D", = Eff\q)jM') - hP(q)Mq') , 
D"; ^ff\q)jt(q')-hP(q)Mq'), 



(3.6) 



with/,(?) = (^+l);>(^)-^;>+i(^)and/f'(^) = {(+\)hP{q)- 
q hpl\iq)- The spherical Bessel and Hankel functions depend on 
q = ka and q' - ■\/sq, with a the radius of the sphere. 

To determine the change in the atomic decay rate due to 
the presence of a dielectric sphere one needs the components 
of the Green function Gc for coinciding arguments. The non- 
vanishing components follow from (13.4b as: 



e^ ■ Gc(r, r,<jj + iO) ■ e,- = 

J CO 

= 7;;i;zT,(-iyv(i+i)fB'^,[hP(kr)f 



(3.7) 



and 



eg ■ Gc(r, r,aj + iO) ■ e^ = e^ ■ Gc(r, r,Lj + iO) ■ e^ 

n2 
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[rh^ikr)] 



dr 



+ B';^[krhf\kr)fy (3.8) 

in agreement with f2l'] . Here e^, eg and e^ are unit vectors in a 
spherical coordinate system. 



4. Decay near a half-space of absorptive scatterers 

We consider a halfspace z < filled with a dilute set of 
spherical scatterers. The non-overlapping spheres are randomly 
distributed with a uniform average density. An excited atom is 
located at r,, = (0,0, Zo), with Za > a, so that the minimal dis- 
tance between the atom and the scatterers is positive. The decay 
rate is given by the sum of the vacuum decay rate and a cor- 
rection term. The modified rate depends on the orientation of 



2 



the dipole-moment transition matrix element. If the dipole mo- 
ment is oriented perpendicular to the z-axis, the vacuum rate is 
Toa = cl>1 Ke|/ij^|§)|^/(37reo^c-')- A similar formula is valid for a 
dipole moment oriented parallel to the z-axis, with fi^^ replaced 

byA^ii- 

If multiple-scattering effects are neglected, the correction 
term in the decay rate is given by the sum of the correction 
terms due to all spheres. For the perpendicular orientation one 
finds 



2a>~ srt 



(4.1) 



•Im GciTa - Ri, Ta - R;, w„ + /O) ■ {g\li^\e) , 

with Gc given by (13.4b and R, the positions of the centers of 
the spheres. Choosing the jc-axis to be parallel to the transition 
matrix element and averaging over the positions of the spheres 
we get 



•^o Jz<0 

X e.r ■ GciVa - r, r„ - r, oja + /O) ■ e.t , 



(4.2) 



with n the uniform density of the spheres and C i a unit vector 
along the .it-axis. The volume integral can be written as a triple 
integral, viz. over |r-rfl|, z and an azimuthal angle. Upon carry- 
ing out the latter two of these integrals one finds for the volume 
integral in (I4.2l i: 





2r2\ 







dr 

Z^ 

~ + -5- ee ■ G,.(r, r, + /O) ■ eg 
3r 3 / 

+ {-Zar + P-) ■ G,.(r, r, + iO) ■ e^] . 
Insertion of ( I3.7l i and (I3.8l l yields 

37mc^ 



(4.3) 



r(uimy(-o'^(^+i) 



4w?, 

x[b?7;^(^„) + b^7»;j^„)] 

with multipole amplitudes given by ( I3.5I )-( I3.6I ) with k 
oja/c, and with the integrals 

^ 2\r.,n,,n2 

l3f3 
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and 



(4.5) 



(4.6) 



with ^ equal to (a - ((^a + iO)Za/c. The derivative of the spheri- 
cal Hankel function in ( 14.5b can be rewritten in terms of Hankel 
functions with a different index ll22ll . For large ( the asymptotic 
forms of these integrals are 



JIJO - (-1) 



2? 



2f 



(4.7) 



so that ( 14.4b becomes 
3nnc^ 



To.^Im^/^ 

(=1 



'+l)[B'c-B'^)e^'^\ (4.8) 



which falls off proportionally to I /(a- Substituting the leading 
terms of Bj, and B™ for small values of q and e - 1 one 
recovers a result found before lllOll . 

Similar expressions may be obtained for the correction to 
the decay rate of an excited atom with a dipole moment parallel 
to the z-axis. Instead of ( 14.21 ) one gets a formula with the zz- 
component of Gc . Upon carrying out the integrals one arrives 
at the analogue of ( 14.41 ). with the integrals 



jy^O = 2e({ + 1) J dt + [hf\t)f 

^r<"(f-f4)[^K"»if 



(4.9) 



and 



dt\^-j^-2^t+—j[h\'\t)] . (4.10) 



For large ^ their asymptotic forms are 

so that the decay rate for large (a becomes 
37Tnc^ 



(4.11) 



(4.12) 



In contrast to ( 14.8b the right-hand side is proportional to the in- 
verse square of f,,. For general values of fa we have to evaluate 
the integrals in ( I4.5b - (l4.6b and ( |4.9b -( l4.10b . as will be done in 
the following section. 



(4.4) 5. Evaluation of integrals 

The integrals J''^ ^ and y^'n are linear combinations of inte- 
grals of the general form 



duu-"h'^^^{^u)h'^^m. 



(5.1) 



which is symmetric in {\,{2- In fact, upon inspecting ( I4.5b -( I4.6I ) 
and ( I4.9b -( l4.10b we find that explicit expressions are needed for 
the integrals I[/^„{0 with n = -2,-1,0, 1,3 and for lee-i.n(0 
for n = -1,0,2. With the use of standard identities ll22ll for 
spherical Hankel functions and by means of a partial integration 
we may derive several relations connecting these integrals for 
different values of the parameters: 



Ie,-l,(2,n(0 + tti + lJiAO - 



2^, + 1 



(0, 



(5.2) 



li + 1 

+(2A + 1) /f„^,+i,„(0 = h^PiO h'l^KO ■ (5.3) 



In order to obtain explicit expressions for with n - 

0, 1, 3 we start from a result [12211 that is valid for n = and 



a, + {2 + i)/f„4,o(f) = 



(5.4) 



as may be checked by differentiation. We omit the argument 
( of the spherical Hankel functions from now on. To obtain 
the corresponding expression for {\ - {2 we put li - i + \, 
I2 - I and n = in ( 15.3b and use (I5.4l i in the second term. 
In this way we obtain a recursion relation connecting /f/ o for 
consecutive values of t. Solving this relation by employing the 
identity /o,o,o(0 ^-(2//0 Ei(-2i() + (with £1 the expo- 



nential integral MM) as an initial condition, we find 
2i 



(2^+ IX 



£i(-2/0 



21 



+ 1 Zj L J 2^ + 1 L M 

k—O 



(5.5) 



for all ( > Q. The exponential integral of purely imaginary ar- 
gument can be expressed in terms of sine and cosine integrals as 
Ei{-2iO = -Ci(2^) - /Si(2^) + in 12. With the help of the iden- 
tities ( 15.2b . (15.4b and the recursion relations for the spherical 
Hankel functions one derives expressions for /f/j (with {£ >\) 
and /f/_3 (with ( > 2), in the form of linear combinations of 
products of spherical Hankel functions: 

r1 



for 



2€{t+\) 
> 1, and 



1 



2^(^-Hl) 2(^-Hl) 



(5.6) 



3(^- 1)^(^ +!)(/' + 2) 6(^ -I- 1)(^ H- 2) 



1 



2{t + 2) 



3i{ -!){({+ l)(£ + 2) 



\){{ + 2) 



2f^ 



\){(+ \){( + 2) 



3(^ + 2) ^^'"^^ 

for t > 2. It turns out that these formulas cannot be used 
for /o,o,i and /i,i,3, since the expressions diverge in these cases. 
However, these special cases can be obtained straightforwardly 
by connecting them to the exponential integral of the same ar- 
gument as in ( 15.5b . 

Expressions for /f/ „ with n - -2 follow by choosing £1 = 
I + \,l2 - i and n - -2 m ( 15.3b . The second term at the left- 
hand side drops out for these values of the parameters. As a re- 
sult a simple recurrence relation for /f/,-2 is found, which may 
be solved for all ^ > by employing the identity /o,o,-2(0 = 
—ie^'^l{2(^) as a starting point. One gets for £ > 0: 

J ()-\ - - - - 21 +\ 



hfh^-^l^. (5.8) 



Furthermore, by choosing in (15.3b the parameters as n = -2 and 
{\ = C2 as either ( or: ( + 1, one gets two identities, which may 
be combined with (15.2b so as to obtain a recursion relation for 
h,c-\- Solving that relation with the initial condition /q _i(^) - 
-C^Ei(-2iO, we find for all ^ > 0: 



/f/,-i(0 = -r'£i(-2/0 + Xi 



jl:=l 



2k +1 
2k(k + 1) 



K'f 



4K"r- 



1 



2(^+1) 



(5.9) 



It should be noted that the sum drops out for f -0. 

Finally, we need expressions for Ie,e-\,p for p - -I and 
p -2. Once more we use the identity (15.3b . now for the choice 
{\ - £2 - ^andn - -1. It yields a recursion relation for 
from which we get for £ >\\ 

t-i 

il:=0 

Turning to the case p = 2, one derives a result for Ie.e-1.2 by a 
repeated use of ( 15.2b in combination with (15.4b . Once again lin- 
ear combinations of products of two spherical Hankel functions 
are found, at least for £ > 2: 



'£,£-1,2 



(0 



3(£- 



!)£{£ + 1) 



3(£-l)£{£+l) 6(£-l) 
2^2 1 

3(£- l)(£+l) 3 



6i£+l) 



"{ "£-l ■ 



(5.11) 



For £ - I this expression is singular and cannot be used. From 
a direct evaluation of the integral for this special case one finds 
that an exponential integral shows up, as before. 

6. Results 

The explicit expressions for the basic integrals (15.1b that we 
derived in the previous section can be employed now to de- 
termine the corrections to the decay rate. For a perpendicular 
orientation of the excited atom the correction (14.4b to the decay 
rate is governed by the integrals (I4.5b - (l4.6b for which we get 
upon substitution of the relevant contributions: 



6£{£+l) 
6£(£+l) 



^\2£^ -2£-3) 

b£{£+l) 
^^(2^2 +2£-3) 



^£ 



2(£ + 1) 



m£ + 1) 

^2(2^2 + 3^ _ 2) 



3(£+ 1) 

£ 



3(^+1) 

3 y _2k+l_ r (,).2 _ 1 . (1). 

^ Zj 2k(k+D L * J 2^ L J 



k=\ 



(6.1) 



and 



6{({+l) 3(^+1) 



2(\2{ + 1) 



3i{+ 1) 



(6.2) 



for > I. Likewise, the results for a parallel atomic orienta- 
tion read 



r,J0^2^Ei(-2iO 



^3(4^2 + 4^ _ 3) 



__2^_ 
1) 

t 



2(\Af + 6t ■ 



+ 1 

2 



1) 



3(^+1) 

n2 



-f'Zl^K'r-f'K'r 



<:=1 



(6.3) 



and 



2^^^(^ - 1) 



3^(^+1) 3(^+1) 

^5 2^3 
'3^(^+ 1) ~ ~ 



2^4 



2^2(2^+1) 



3(^+1) 



j(:=l 



(6.4) 



again for all ^ > 1 . As remarked above, the exponential inte- 
grals can be expressed in sine and cosine integrals. 

After insertion of the expressions (I6.H - (I6.2) and the mul- 
tipole amplitudes ( 13.5b into ( I4.4l i. the average correction to the 
decay rate for the perpendicular configuration is found in terms 
of well-known special functions depending on („, q and s. It 
may be plotted as a function of (a, for various choices of q 
and e. To facilitate comparison with our previous results llioll 
we introduce the decay rate correction function f±i(a, q, s) = 
-16(rc,i)/(3nvoro,i) with vq - Ana^ /3 the volume of the spheres. 

We shall first consider two special cases that we have treated 
in iflOil : purely scattering spheres and purely absorbing spheres. 
In the former case we choose the dielectric constant to be real 
{s - 1 .5) and the spherical radius to be finite on the scale of 
the wavelength (q = 0.5). For the purely absorbing case with 
vanishingly small spheres (q — > 0), we take the dielectric con- 
stant to be complex with the value e = 1.5 H- /0.5, as in |10]. 



Since for small q the multipole amplitudes B'j. and B'" behave as 
qie+i gr2^+3^ respectively, only the electric dipole amplitude 
Bj = iq^(s- 1)/(e+2) contributes to (14.41 1 for the purely absorb- 
ing case. In Figs. [T] and |2] the curves for f±{^a) are compared 
to their asymptotic counterparts for large (a that follow from 
(14.81 1. As can be seen from these figures, the asymptotic curves 
are quite adequate already for ~ 3. In the asymptotic regime 
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Figure 1: Decay rate coiTection function fi_{{a) (solid line) and its asymptotic 
form at large distances (dashed line), for a medium with scattering spheres (with 
q = 0.5,£(cja)= 1.5). 
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Figure 2: Decay rate correction function f±((a) (solid line) and its asymptotic 
forms at small and large distances (dashed lines), for a medium with absorbing 
spheres (with q = 0, e(Wo) = 1.5 + (0.5). 



the results given in II Oil are corroborated. (It should be noted 
that the curves given in ifioll erroneously represent - fj_ instead 
of ./±.) For small distances the behaviour of the atomic decay 
rates in the two cases differ considerably. In fact, for the purely 
scattering case of Fig. [T] the decay rate attains a finite value 
when approaches its minimum value q. On the other hand, 
for the purely absorbing case of Fig. |2]the decay rate correc- 
tion function is governed by J''^ ^{0, which according to (16. Il l 
has the asymptotic form -1/(4^^) for small ^. Hence, the decay 
rate correction function diverges as -(3 / 2) Im[(£- 1 ) /(e+2)] /^^ 
for — » in this case. 

For a more general situation in which both scattering and 
absorption take place, we choose q = 0.5 and s - 1.5 + /0.5, 
with results presented in Fig.|3] For large (a the decay rate falls 
off" like in agreement with ( I4.8l l. For small distances the 
decay rate diverges, as in Fig. |2] In fact, as (a — > q, the rate 
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Figure 3: Decay rate correction function f±{^a) (solid line) and its asymptotic 
forms at small and large distances (dashed lines), for a medium with scattering 
and absorbing spheres (with q = 0.5, £(<jja) = 1-5 + (0.5). 



Figure 4: Decay rate correction function /i|(fa) (solid line) and its asymptotic 
forms at small and large distances (dashed lines), for a medium with scattering 
and absorbing spheres (with q = 0.5, e(<jJa) = 1.5 + /0.5). 



is found to be proportional to ((a - ^) This follows from 
( 14.4b . since the series converges increasingly slowly when (a 
approaches q. Indeed, for large i the electric multipole ampli- 
tudes are given by [i^/(l^[(2{- l)!!]^)] [(g- l)/(e+ 1)] q^^+\ 
while the integral dO gets the form -[(2{ - l)!!]2/(2^f 
Hence, the electric multipole contribution to the {-th term in the 
series of i4~^ is - l)/(e + 1)] (ql(,af^^\ Since the mag- 

netic multipole contributions turn out to be negligible for large 
{, the asymptotic form of ( 14.41 1 for ^„ tending to q is propor- 
tional to (^/^fl)^'^^' - ql{2{^„ - q)}, so that the asymptotic 
form of /j. reads 



Im 



s- 1 



s+ 1 



(6.5) 



for (a q. 

The physical mechanism for the divergence in the decay 
rates of Figs. |2]and|3]for small is the efficient non-radiative 
energy transfer from the atom to the absorbing spheres that 
dominates the atomic decay in the near zone. A similar diver- 
gent behaviour has been found in a classical treatment of the en- 
ergy transfer between an excited molecule and a homogeneous 
absorbing medium filling a halfspace [23]. It should be noted 
that ( 16.51 ) loses its meaning when (a - q becomes so small that 
the approximations made in deriving it are no longer valid. In 
particular, perturbation theory in lowest order and the electric- 
dipole approximation are not adequate to describe the decay for 
very small distances. Furthermore, the notion of scatterers with 
a structureless surface gets lost as well in that case. 

For the parallel configuration we have likewise evaluated 
f\\((a,q,s) = -16(rc,||)/(3nvor().||). The result for the mixed 
case of both scattering and absorption is given in Fig.|4]for the 
same choice of the parameters q and s as in Fig.|3] The figure 
clearly shows that for the parallel configuration the correction 
to the atomic decay rate goes faster to zero with increasing 
than for the perpendicular configuration. This is in accordance 
with the findings of section 4, where it has been seen that in 
the asymptotic regime the correction to the atomic decay rate 
is proportional to the inverse distance in the perpendicular con- 
figuration, but to its square in the parallel configuration. As 
before, the asymptotic expression is adequate from ~ 3 on- 



wards. For small distances, with (a — » q, the asymptotic form 
of /ii is twice that of f±, as follows by comparing the asymptotic 
forms of ( 16. Il l and (I6.3l l for large d. 

In conclusion, we have shown how absorption and scatter- 
ing processes in a medium may cooperate in modifying the 
emission rate of an excited atom in its vicinity. The explicit 
expressions for the decay rate that we have obtained permit a 
detailed analysis of the behaviour of the emission rate for arbi- 
trary distances between the atom and the medium. As we have 
seen, the effects of absorption and of scattering are qualitatively 
different, when the atom approaches the medium. 
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